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Reconstructing Polyatomic Structures from Discrete X-Rays: 
NP-Completeness Proof for Three Atoms 

Marek Chrobak* Christoph Diirr^ 



Abstract 



f~| ' We address a discrete tomography problem that arises in the study of the atomic structure 

^ , of crystal lattices. A polyatomic structure T can be defined as an integer lattice in dimension 

D >2, whose points may be occupied by c distinct types of atoms. To "analyze" T, we conduct 

i measurements that we call discrete X-rays. A discrete X-ray in direction ^ determines the 

^^ ' number of atoms of each type on each line parallel to ^. Given £ such non-parallel X-rays, we 

wish to reconstruct T. 
^Q ' The complexity of the problem for c = 1 (one atom type) has been completely determined 

^^ , by Gardner, Gritzmann and Prangenberg |^, who proved that the problem is NP-complete for 

any dimension D > 2 and ^ > 3 non-parallel X-rays, and that it can be solved in polynomial 
QJ i time otherwise Q. 

The NP-completeness result above clearly extends to any c > 2, and therefore when studying 
the polyatomic case we can assume that i = 2. As shown in another article by the same authors, 
^ ■ W, this problem is also NP-complete for c > 6 atoms, even for dimension D = 2 and axis-parallel 

OO ! X-rays. The authors of M conjecture that the problem remains NP-complete for c = 3,4,5, 

although, as they point out, the proof idea in H] does not seem to extend to c < 5. 

We resolve the conjecture from [El by proving that the problem is indeed NP-complete for 
c > 3 in 2D, even for axis-parallel X-rays. Our construction relies heavily on some structure 
^^ \ results for the realizations of 0-1 matrices with given row and column sums. 

Y^ ' Keywords: discrete tomography, high-resolution transmission electron microscope, multi-commodity 

• • . max flow. 
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1 Introduction 

The fundamental principle of the transmission electron microscope (TEM) is very similar to the 
more familiar optical microscope: it "shines" a focused beam of electrons towards a specimen, 
and the transmitted beam is projected onto a phosphor screen, thereby generating an image. The 
intensity represents the density and thickness of the specimen: denser or thicker areas of the 
specimen transmit fewer electrons and produce darker areas in the image. The development of the 
TEM in 1930's was necessitated by the limitations of the optical microscopes, whose magnification 
and resolution were insufficient to study the internal structure of organic cells or to find defects 
in bulk materials. Recently, new advancements in high-resolution TEM (HRTEM) led to the 
development of instruments and techniques for studying biological molecules and for investigating 
the atomic structure of crystals. In particular, a technique called QUANTITEM [Q, IC] allows 



one to determine the number of atoms in the atom columns of a crystal in certain directions. 
Given these numbers, we wish to reconstruct the structure of the crystal. This is an example of 
an algorithmic problem belonging to discrete tomography, the area of mathematics and computer 
science that deals with inverse problems of reconstructing discrete density functions from a finite 
set of projections. The size of crystals that occur in materials science applications is about 10^ 
atoms, and, for data sets that large, efficient reconstruction algorithms would be of great interest. 

The problem we address in this paper can be formulated as follows: Define a polyatomic structure 
T as an integer lattice in dimension D > 2, whose cells may be occupied by c distinct types of 
atoms. Each of these cells can be occupied by one atom, or it could be empty. To "analyze" T, we 
conduct i measurements that we refer to as discrete X-rays. (QUANTITEM uses electron beams, 
but, following ||5[, we use a more familiar term "X-ray" instead.) A discrete X-ray in direction ^ 
determines the number of atoms of each type on each line parallel to ^. Given such i non-parallel 
X-rays, we wish to reconstruct T. 

The complexity of the problem for c = 1 (one atom type) has been completely determined 
by Gardner, Gritzmann and Prangenberg Q, who proved that the problem is NP-hard for any 
dimension D > 2 and i > 3 non-parallel X-rays, and that it can be solved in polynomial time 
otherwise [^. 

The NP-hardness result above clearly extends to any c > 2, and therefore when studying the 
polyatomic case we can assume that £ = 2. As shown in another article by the same authors, Q, 
this problem is also NP-hard for c > 6 atoms, even for dimension D = 2 and for the axis-parallel 
X-rays. The authors of Q conjectured that the problem remains NP-hard for c = 3, 4, 5, and they 
pointed out that for these values of c "a substantially new technique will be needed, at least for 
the case c = 3" . 

We resolve the conjecture from Q by proving that the problem is indeed NP-hard for c = 3 
(and thus for any larger c as well) in 2D, even for the orthogonal case, that is, with axis-parallel 
X-rays. 

In the orthogonal case, the problem is equivalent to that of reconstructing (c-l-l)-valued matrices 
(c atom types and "holes" ) from the row and column sums for each atom. Without loss of generality, 
we can concentrate on square, say L x L, matrices. Let A be the set of c atom types. For any atom 
type a G A, denote by rf (resp. s^ ) the row-sum (resp. column-sum) of atom a, that is, the number 
of atoms of type a in row i (resp. in column j). The vectors r" = (rf, . . . , r£) and s'* = (sf , . . . , s^) 
are referred to, respectively, as the row-sum vector and the column-sum vector for atom a. 

A realization of the sums Z = (r°, s")agA is an L x L matrix T with values from A U {□}, such 



that for each atom type a G A 
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We say that X is consistent if it has a reahzation. 

More specificahy, we concentrate on the following decision problem: 

c-Color Consistency Problem (c-CCP) 

Instance: row and column sums I = (r"',s"')a(zA, where |A| = c; 
Query: Is I consistent? 

Gardner, Gritzmann and Prangenberg proved in Q that 6-CCP is NP-complete. In this paper 
we prove that 3-CCP is NP-complete. 

If we restrict ourselves further to just one atom (that is, 1-CCP), the problem becomes equivalent 
to the reconstruction of 0-1 matrices from the row and column sums - a problem predating the 
discrete tomography research. The first efficient reconstruction algorithm was proposed in 1963 
by Ryser Q, and a similar algorithm was rediscovered in 1971 by Chang [||. In addition to 
reconstruction, Ryser and others studied various structural properties of 0-1 matrices with given 
row and column sums, and our construction relies heavily on some results in this area. Interested 
readers are referred to an excellent survey by Brualdi |l[. 

The general idea of the proof is explained in Section ^. In Section ^, we review the structural 
properties of 0-1 matrices with given row and column sums that are needed for our proof. Using 
these properties, we prove the Skew-Mirror Lemma in Section ^. In Section |5|, we construct a 
number of gadgets, including "skew mirrors" and "edge verifiers" , and we prove that they satisfy 
the desired properties. Finally, in Section y, we present the complete construction and give the 
formal NP-completeness proof. 

In addition to the QUANTITEM method, the problem of reconstructing lattice sets from their 
projections arises naturally in a variety of other areas, including statistics, data security, and image 
processing. It can also be expressed as a multicommodity flow problem. We discuss these issues in 
Section |^ where we also comment on the last unresolved case, c = 2. 

2 The General Idea of the Proof 

In the proof, we use a reduction from the Vertex Cover problem: 

Vertex Cover Problem 

Instance: An undirected graph G(y,E), an integer K; 
Query: Is there a vertex cover of G of size K? 

Recall that a vertex cover of a graph G = {V, E) is a set U C.V such that for all (u, v) S E, 
either u S C/ or v E [/. The Vertex Cover problem is well known to be NP-complete (see, for 
example, ^). Let n = \V\ and m = \E\. We assume, without loss of generality, that m,n> 1. 

Suppose first that, using some set A' of d atom types, we can force a unique realization of the 
form shown in Figure ||. We call this realization a frame. In the frame, the empty entries form two 
diagonals, the main diagonal of length (m + l)n, and the side diagonal of length mn. We divide 
both diagonals into intervals of length n that we refer to as mirrors. Thus we have two rows of 
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mirrors: m + \ mirrors in the main-diagonal, and m mirrors in the side-diagonal. All other entries 
are filled with atoms from A'. 

We now add two more atom types C,D ^ A'. Use 
atom D to create m copies of a candidate vertex cover 
U in the following way: The first row and column D- 
sum is K and all other D-smns are 1. (Figure || shows 
the D-sums.) Then the pattern of Ds in each main- 
diagonal mirror is the same, and is also the same as 
the pattern of holes in the side-diagonal mirrors. We 
associate U with this pattern: a vertex u is in [/ iff 
the uth cell in any side-diagonal mirror is a hole. We 
think of [/ as a "beam" projected onto the last n cells 
in the first column, repeatedly reflected in a double-row 
of mirrors, and exiting through the last n cells in the 
first row. 

Finally, we can use atom C to verify that U is indeed 
a vertex cover. In order to do so, we convert the jth 
side-diagonal mirror into an edge verifier for edge ej = Figure l: The frame and mirrors for m^ 3. 
(n, v) (it may be necessary to add some more rows and columns to the matrix shown in Figure IT]). 
Using appropriate sums for atom C, the realization of atoms in A' can be extended to a realization 
of all atoms, including C, iff either the uth cell or the uth cell in side-diagonal mirrors is a hole 
(and thus, either n E C/ or f G [/). 

An idea similar to the one described above was used by Gardner, Gritzmann and Prangenberg 
Q (they used a reduction from a different problem, not Vertex Cover). Using 4 atoms they 
constructed, in essence, what we call a frame, obtaining the NP-completeness proof for 6 atoms. 
In our first attempt to improve their construction we were able to construct the frame gadget with 
only 3 atoms, reducing the total number of atoms to 5. However, this idea does not work when 
fewer than 5 atoms are available. As pointed out by Q, a new approach is needed. 

The main idea behind our proof is this: Define a partial order "^" on all iC-element vertex sets 
(candidate vertex covers). The important property of "^" is that its depth is polynomial, namely 
at most J = K{n — K) + 1 (each strictly increasing chain has length at most J). Further, "^" has 
a unique minimum element [/""", and a unique maximum element f/"""^. Instead of using "perfect" 
mirrors, we use "skew" mirrors. These mirrors have the property that the reflected set is never 
smaller (with respect to partial order "^") than the set projected onto a skew mirror. These skew 
mirrors are also "wobbly" — we know that they can reflect the same or a bigger set, but we cannot 
control what exactly the reflected set will be. 

Now, instead of using m mirrors, we use mJ skew mirrors in the side-diagonal. They are divided 
into J segments of m mirrors each. In each segment, the jth skew mirror is converted into an edge 
verifier for edge Cj. We "shine" [/'"'" onto the first mirror in the bottom-left corner, and we make 
sure that the final set resulting from all reflections in the top-right corner is f/'""^. Since "^" has 
depth J, there has to be a segment in which all mirrors reflect the same set U . Then the edge 
verifiers in this segment will verify that U is indeed a vertex cover. 

Why does it help? It turns out that our skew mirrors can be constructed using only two 
atom types. Furthermore, the same atom types can be used to encode the information about the 
candidate vertex cover U . We use one more atom type to construct edge verifiers, and thus we only 
need three atom types for the whole construction. 



3 0-1 Matrices with Given Row and Column Sums 

In this section we review some basic results from the Hterature on 0-1 matrices with given row and 
column sums. 

By X, y, z we denote nonnegative integer vectors of length p, for example x = (xi, . . . , Xp). The 
reconstruction problem for 0-1 matrices with given row and column sums is equivalent to 1-CCP, 
and can be stated as follows: Given x and y, is there a 0-1 matrix T that has Xi I's in row i and 
Uj I's in column j, for all 1 < i,j < p? Again, in this case, a matrix T satisfying these conditions 
is called a realization, and x, y are called consistent if they have a realization. 

The structure function. Given a. p x p matrix T, and integers < k,l < p we partition T 
into four submatrices (which may have zero width or height): Tj^, T^, T^ and T^^, defined by the 
intersections of the first k rows or the last p — k rows with the first / columns or the last p — I 
columns), that is 



T 






By |T|^ and \T\q we denote the numbers of I's and O's in matrix T. 
For a given instance x, y, the structure function t/^i is defined by 



/ p 

Tki = {p - f^){p - ^) +^yj - J2 ^■'■ 



j=l i=k+l 



Then for any arbitrary realization T we have 



I p 

Tki = {p - k){p - 1) +^yj - J2 ^* 

j=l i=k+l 
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Consistent sums. We now show that, using the structure function, it is possible to characterize 
consistent sums. An integer vector z = (zi, . . . , Zp) will be called monotone if zi < . . . < Zp. 

Lemma 1 Q Monotone vectors x,y are consistent if and only if t^i > for all k,l = 1, . . . ,p. 

The implication (=>) in Lemma || follows directly from Equation (||). The implication (<^) can 
be proven constructively by giving an algorithm that produces a realization T for any pair (x, y) 
for which the structure function is non-negative. (See ||l| for details.) It is also not hard to see that 
Lemma |^ can be derived from the Max-Flow-Min-Cut theorem for network flows. 

Decomposed realizations. We say that T is (k, I) -decomposed if T^ consists only of O's and 
T^ consists only of I's. The lemma below follows immediately from Equation (|^), and it will play 
a major role in this paper. Note that in this lemma, as in the definition of r, we do not require the 
projections x,y to be monotone. 

Lemma 2 Q Suppose that T is a realization o/x, y, and let < k,l < p. Then Tki = if and 
only ifT is (k,l)-decomposed. 



Remark 1 Lemma^ implies that if just one realization o/x, y is {k, I) -decomposed, then all real- 
izations are {k, I) -decomposed as well. 

4 The Skew-Mirror Lemma 

0-1 Vectors and minorization. We use Greek letters a,(3, . . . for 0-1 vectors of length p, say 
a = (tti, . . . , Op). The complement a of a is Sj = 1 — ckj for all i = 1, . . . ,p, and the reverse a is 
5i = ap-i+i ioT i = l,... ,p. 

We say that a minorizes /3, denoted ct ^ /3, if 

k k 

^ai < ^Pi yk = l,...,p. 

i=l i=l 

By straightforward verification, "^" is a partial order. We also write a ^ (3 if ex ^ P and a ^ fi. 
The tota/ sum of a 0-1 vector a is Eq = X]f=i c^j- If '^j /^ a^^e two 0-1 vectors with equal total 
sums, then the definitions above imply directly the following equivalences: 

cx^P ^=^ a ^(3 ^^ ahP- 

An important property of the minorization relation is that it is "shallow", that is its depth 
is only polynomial (unlike, for example, the lexicographic order). The next lemma gives a more 
accurate estimate on the depth of "^". 

Lemma 3 Suppose that we have a strictly increasing sequence of 0-1 vectors 

Oi^ -< o? ^ . . . ^ oP, 

with total sums Ea* = t for each i. Then q <t{p — t) + 1. 

Proof To each 0-1 vector a. assign the number ||q;|| defined by 

p p k 

j=l k=l i=l 

If a ^ /3 then Yli=i ^i — Si=i ft for A; = 1, ... ,p, and this inequality must be strict for at least 
one k. We conclude that a ^ f3 implies ||q;|| < \\(3\\. 

Now, by the argument above, the numbers ||q!*|| are strictly increasing. Therefore 

q < IIq^II - ila^ll -Fl 
p t 

< E k-^k + l 

k=p-t+l k=l 

= t{p-t) + l, 

completing the proof. ■ 
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(a) (b) 

Figure 2: (a) A realization of x, y for a ~ 010100 and /3 — 000101. Disks represent I's. The row and 
column sums that correspond to O's in a and (3 are shown in hold, (b) Perfect mirror PMa- for cr = 101011. 
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The 0-1 skew mirror. The lemma below deals with a special instance of the reconstruction 
problem for 0-1 matrices, in which the row sum vector x is determined by a 0-1 vector a, and the 
column-sum vector y is determined by a 0-1 vector f3. 

Given a 0-1 vector a of length p, we associate with cr apxp perfect mirror matrix PM^- defined 

by 



PM^[i,j] 



In a perfect-mirror matrix the cells on the main diagonal i + j = p + 1, counted from top down, 
contain a, while all cells above it are 0, and all cells below it are 1 (see Figure |2|b). From Lemma y 
we immediately obtain the following corollary. 

Corollary 1 Let a be a 0-1 vector of length p. Then PM^- is a realization of vectors x, y if and 
only if for each k = 1,... ,p, (a) Tk,p-k = 0, and (b) Tk,p+i-k = iff ak = 0. 

Note that Corollary || together with Remark |l| implies that if PM,,- is a realization of x, y then 
it is unique. 

Lemma 4 Let a,P be two 0-1 vectors of length p, and letx,y be row and column sums defined by 
Xi = i — Oi and yi = i — Pi, for i = 1, . . . ,p. Then 

(a) Vectors x, y are consistent iffT,a = Ti(3 and a^ (3. 

(b) Suppose that x,y are consistent. Then a = f3 iff the unique realization o/x, y is PM^- 

Lemma ^ is illustrated in Figure Q. In Figure ^(a), we have a = 010100, and (3 = 000101 
(we write 0-1 vectors as binary strings, for simplicity), for which the corresponding row sums 
and columns sums are x = (1,1,3,3,5,6) and y = (1,2,3,3,5,5). Note that S5 = S/3 = 3, 
and that 5^/3, implying that x and y are consistent. One realization of x, y is shown in 
Figure |2|(a) . In Figure 0(b), a and x are the same as in (a). Since y = (1,2,2,4,4,6) corresponds 
to /3 = 001010 = 5, vectors x, y have a unique realization that is a perfect mirror PM^-, for 
a = 101011 = a. 



Proof Vectors x, y are monotone, so we can use Lemma ||. We start by computing the structure 
function for x and y: 

I p 

Tkl = {P - k){p - I) +J2yj ~ ^ ^i 

jr=l i=k+l 

I p 

j=l i=k+l 

P I 

= E «.-5]/?,- + i[2(j,-A;)(p-/)-p(p + l) + A:(A; + l) + /(/ + l)] 

i=k+l j=l 

p—k I 

= ^^i-^Pj + Up-k-l-l){p-k-l). (2) 

Now we are ready to prove Part (a). We prove the two imphcations separately. 

(=>) For any / = 0, ...,p, using Equation (||) with k = p — I, we get that Tp_i^i > imphes 
^i=i Oii — Yl,j=i l^j — 0- Thus a ^ 13. Moreover, x and y have equal total sums, if and only 
if Eq = S/3. 

(■<=) Assume that 5^/3. We consider two cases, when k + I < p and k + I > p + 1. 
Suppose first that k + I < p. From Equation (^) we have 

I p-k I 

Tkl = y^5i+ y^ ai-y^^f3j + ^{p- k-l -l){p- k-l) 

i=l i=l+l j=l 

> 0, 

because Yli=i ^« ~ Y2j=i l^i — 0' ^^^ {p — k — I — l){p — k — I) > 0. 
Suppose now that k + I > p + l- From Equation (^) we have 

I I I 

Tkl = E "* ~ E o^i-'^Pj + ^{p-k-l-l){p-k-l) 

i=l i=p—k+l j=l 

I I I 

= E"i~E^J~ E Oli + ^{k + l + l-p){k + l~p) 
i=l j=l i=p—k+l 

I I I 

= E"^~E'^i + ^ E {k + l + l-p-2ai) 

i=l j=l j=p— fc+1 

> 0, 

because X]j=i ^j ~ ^12^=1 Pj — 0' ^^^ k + l + 1— p — 25j > 2 — 2Sj > 0. 

Now we prove Part (b). By Corollary |l| and Equation (pi), a realization of x, y is a perfect 
mirror PMa-, for some cr, if and only if 5 = /3. Thus it is sufficient to show that a = (3 implies 
that (7 = a. This follows by simple verification of row sums. ■ 



5 Some Useful Gadgets 

Recall that G,K is the given instance of Vertex Cover, where G = {V,E), n = \V\ and m = \E\. 
Without loss of generality we can assume that < K < n. All our examples will refer to the graph 
of Figure y. In this figure, G has n = 6 vertices, m = 3 edges, and set {1,3, 5, 6} is a vertex cover 
with K = 4 vertices. (Clearly, this is not a minimum-size vertex cover. We use this example only 
to illustrate certain features of our construction.) 



Figure 3: Example of a vertex cover of size K = 4. 

The purpose of this section is to introduce two (n + 2) x (n + 2)-instances of 3-CCP called the 
beige skew mirror and the edge verifier. They will be used later in the NP-completeness proof. 

Throughout the rest of the paper we will use capital letters A, B, C to denote the three atom 
types, and we will sometimes refer to these types as colors: Azure, Beige, and Cyan. 

Beige skew mirror. Given two 0-1 vectors cx,(3 of length n, we define the beige skew mirror as 
an (n + 2) x (n + 2) instance of 3-CCP, BSM(a, (3) = (x^, y^), with the following row and column 
sums: 

= i — Pi + 2 for i = 1, . . . , n 

= n + 2 for i = n + 1,71 + 2. 

The azure and cyan sums are zero. Figure Q(a) shows an example of a beige skew mirror. 

Lemma 5 Let a, (3 be two 0-1 vectors of length n. Then BSM{a,P) is consistent if and only if 
T,a = Tj(3 and a y (3. 

Proof By definition, any realization of BSM{a, (3) has its last 2 rows and last 2 columns com- 
pletely filled with beige atoms. Define 
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11 



where i = 1, . . . ,n. Then BSM{a, f3) is consistent if and only if the instance (x, y) of 1-CCP is 
consistent. Applying Lemma Q, we obtain that BSM[a, P) is consistent if and only if Sec = S/3 
and a y f3. m 

Azure skew mirror. Given two 0-1 vectors 7, S of length n, we define the azure mirror as an 
(n + 2) X (n + 2) instance of 3-CCP, ASM{'y,6) = (x^,y^, x^,y^), with the following row and 
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A = 0, 
B = O and 



Figure 4: (a) A realization o/ B5'M(010100, 000011). (b) A realization o/ SF(010100, 001010, (3, 5)); it 
verifies that the vertex set {1,3,5,6} covers the edge (3,5). Solid filled regions represent entries that are 
filled with beige atoms independently of the parameters of BSM and EV. 
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for i = 1, 



for i = 1, 
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Lemma 6 Let 7, S be two 0-1 vectors of length n such that S7 
is consistent if and only if ^ ^ S . 



E(5 



n 



K. Then ASM{j,S) 



Proof We claim that each realization of ASM{-f, 5) has beige atoms on positions: 



{ii,n + 2):ji = l} 
{(n + 2,i) ■.5i = l} 
{(n + l,n+l),(n + l,n + 2),(n + 2,n + l),(n + 2,n + 2)}. 



(last column) 

(last row) 

(lower right 2x2 corner) 



Since x^^2 = n — K + 2, and there are exactly n — K -\- 2 non-zero beige column sums, the beige 
atoms are determined uniquely, as shown above. Similarly, the beige atoms in the last column are 
determined uniquely. The last yet unallocated beige atom must be at (n + l,n + 1). 

We now examine azure atoms. Row n + 1 and column n + 1 have no azure atoms. In row n + 2 
azure atoms are forced to be in columns i for which 6i = 0, since all other positions are occupied 
by beige atoms. Similarly, in column n + 2 azure atoms are in rows i for which 7^ = 0. 

Let X, y be the following row and column sum vectors: 



Xi 

Vi 



t-li 
i - 6i, 



for i = 1, . . . ,n. Then ^5'M(7,5) is consistent iff the instance (x,y) of 1-CCP is consistent. By 
Lemma Q, this is equivalent to 7 ^ 3, or 7 ^ 6, completing the proof. ■ 
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Edge verifier. For two 0-1 vectors 7, S of length n, and for an edge e = (u, v) (with u < v) we 
define the edge verifier for e, as a (n + 2) x (n + 2) instance of 3-CCP, 



EVij,S,e) = (x^y^x^y^x^,y^), 

where the azure and beige sums are exactly the same as in the azure skew mirror ASM(-y, S), and 
the cyan sums are: 



< - 
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= 1. 



Lemma 7 Let 'y be a 0-1 vector of length n, and e = {u, v) (with u < v) be an edge of G. Then 
EV{-f, 7, e) is consistent if and only if E7 = n — K and either 'ju = or 7^ = 0. 

Lemma ^ has the following interpretation: if we associate with 7 the vertex set U = {u : ju = 0}, 
then EV{'-f,^,e) is consistent if and only if at least one endpoint of edge e belongs to U. See 
Figure Q(b) for an example of an edge verifier. 

Proof We can assume that S7 = n — K. By Lemma ^, ASM{'y, 7) is consistent. Furthermore, 
by Part (b) of Lemma ^, for any 1 < i,j < n, a realization T of ^S'M(7, 7), satisfies: T'[i,j] = □ 
for i + j < n, T[i, j] = A, for i + j > n + 2, and for i-\- j = n + 1 we have the following equivalence: 
r[i,j] = niff7, = 0. 

If EV{-f,^,e) is consistent, we can extend T to a realization of EV{-f,^,e), and consider the 
positions of cyan atoms. The position (n + l,n + 1) contains a beige atom and {v,n — u + 1) an 
azure atom. This leaves these possible positions for the cyan atoms: 

{u,n — V + 1) (n, n — u+1) (u,n+l) 

{v,n — v+l) (f,n+l) 

{n-\- l,n — V + 1) {n-\- l,n — u + 1). 

We claim that either T[u, n — u + 1] = C or T[v, n — v + 1] = C. For otherwise, the cyan row 
sums x^ = 2 and Xy = 1 force T[u, n + 1] = T[v, n + 1] = C, contradicting "Un+i = 1- 

In summary, we get that £"^(7,7,6) is consistent iff one of 7^, 7^, equals 0. ■ 



6 The Proof of NP-Completeness 

In this section we give the overall reduction. We will define an instance of 3-CCP, which has a 
solution if and only if there is a vertex cover. Its matrix can be divided into regular quadratic 
blocks and we will show that in every realization some blocks will be realizations of particular beige 
skew mirrors, or of particular edge verifiers, while the remaining blocks are filled either with beige 
or azure atoms. For this purpose we partition the matrix into regions of different shapes, which we 
call mirrors, gutters, screens and frames, and we show some particular properties of them. Then 
we conclude that these properties imply the existence of a vertex cover. 
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6.1 The Reduction 

Recall that G = (y,E),K is the given instance of Vertex Cover, where \V\ = n, \E\ = m and 
< K < n. Define J = K{n — K) + 1 and L = (mJ + l)(n + 2). We now show how to map G, K 
into an L X L instance of 3-CCP 

I = (r^s^r^s^r^,s^). 

To specify the sums in Z, it is convenient to view L x L-matrices as being partitioned into (mJ+ 1)^ 
submatrices of size (n + 2) x (n + 2), called blocks. A row or column index is then defined by its 
block index a = 0, . . . , mJ and offset i = 1, . . . , n + 2. For a ^ mJ the azure and beige sums are: 











i 


i = 1, . . . ,n 


' a{n+2)+i 


- „A 

- *a(n+2)+J - 


= (mJ- a- l)(n + 2) + 
r i + 2 


< 



K 


i = n + 1 
i = n + 2 

i = 1, . . . ,n 


' a{n+2)+i 


- *a(n+2)+J - 


= a(n + 2)+^ n + 4 

1 2n + 4- 


K 




i = n+1 
i = n + 2 



and for a = mJ the azure sums are zero and the beige sums are 



i + 2 i = l,...,K 

''a{n+2)+i ~ *'a(n+2)+i ~ " " 



ra(n+2)+i = 4(n+2)+i = a{n + 2) + { i + 1 i = K + l,...,n 



n+2 i = n+l,n + 2. 

Finally, we define the cyan sums. For j = 0, . . . ,J — 1 and k = 0, . . . ,m — 1 let a = jm + k and 
b = mJ — 1 — a. If Ck = (n, v) (with u < v), then 



a(n+2)+u 


= 2 


''6(n+2)+n-M+l 


= 1 


a{n+2)+v 


= 1 


'^b{n+2)+n-v+l 


= 2 


a{n+2)+n+l 


= 1 


„C 
*b(n+2)+n+l 


= 1 



The row and column sums not defined above are assumed to be zero. 

6.2 Realizations of Azure and Beige Atoms 

Let A and S be (n + 2) x (n + 2) matrices completely filled with azure and beige atoms, respec- 
tively. We will use notation A{'y,S) for realizations of ASM{-f,S) and B{a,P) for realizations of 
BSM{a,(3). We define tv by 

TT = . . . 1 . . . 1 . 

K n-K 

For 0-1 vectors P , a^, . . . , P"^ ~ , ct™ , each of total sum n — K, consider L x L azure-and-beige 
matrices of the following form: 
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A 


A 




A 


A 


A 




A 


A 


A 




A 


A 


A 




A{a^f3') 


A 


A{cy\ 


0') 


B{oc\(3^) 


4K,/3") 


B{cP, 


a') 


B 


B(7r,/3") 


B 




B 



B B 



B 
B 
B 
B 



B 
B 
B 

B 



(3) 



Lemma 8 Let I^^ he the restriction of X to the azure and beige sums only. Then a matrix T is a 
realization of Z^^ if and only if T has the form ^, where 



<^o 



^1 



7r^/3 ^a" ^13 < cx^ < . . . ^ (3 



mj—l 



< a 



^ TT. 



(4) 



Proof Note that by Lemma |5| and y, a matrix of the form (y) exists iff inequahties (Q) are true. 
(<^) Let T be a matrix of the form (P). By straightforward verification of the row and column 
sums we obtain that T is a realization of X^^ . (Note that the entries of Q;*,/3*, i = 0, . . . , mJ — 1 
appear in the beige sums with the plus sign in the azure skew mirrors, and with the minus sign in 
the beige skew mirrors.) 

(=^) Let now T be a realization of T"*^, and denote by F an arbitrary matrix of the form (^). Block 
(a, 5) consist of entries in rows a{n + 2) + i and columns h{n + 2) + j for all i, j = 1, . . . , n + 2. We 
call it 



an upper-left block 
the side-diagonal block b 
the main- diagonal block b 
and a lower-right block 



if a + 6 < mJ, 
if a + 6 = mJ, 
if a + 6 = mJ + 1, 
if a + 6 > mJ + 1. 



We claim that T has the structure depicted in Figure |5[ Let T^ be the 0-1 matrix representing 
the non-azure cells in T: T^[i,j] = 1 if and only if T[i,j] ^ A. Let F^ be the analogous matrix 
for F . Then T^ and F^ are realizations of the same instance of 1-CCP. Therefore, if F^ is (A;, /)- 
decomposed so must be T^ . Let k = [mJ — d)(n + 2) and / = d{n + 2) for any d = 0, . . . , mJ. 
Then F^ is {k, /)-decomposed. Thus all upper-left blocks in T are exactly A and all main-diagonal 
and lower-right blocks have no azure atoms. The union of all upper-left blocks will be called the 
A-frame. 

Let T^ be the 0-1 matrix T representing the beige atoms in T, that is T^[i,j] = 1 if and only 
if T[i,j] = B. Let F^ be the analogous matrix for F. Let k,l be a pair of indices such that 
either k = {mJ -|- 1 — d)(n + 2) — 2 and / = d{n + 2) — 2 for some d = 1, . . . , mJ, or (k, I) G 
{(0, L — 2), (L — 2,0)}. Then, by the structure of realizations of azure and beige skew mirrors in 
Lemmas g and g, F^ is {k, /)-decomposed, and so must be T^. The region in T corresponding to 
the I's in the submatrices Tj^ (for the k,l chosen above) is called the B-frame. So the B-frame of 
T is all beige. Since the lower-right blocks are included in this region, they are exactly B. 
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A-frame 



A-mirror 



h-gutter 
h-screen 



B-miiTor 



A-mirror 1 



h-gutter 1 
h-screen 1 



B-mirror 1 



A-mirror 2 



h-gutter 2 



h-screen 2 



B-mirror 2 



B-mirror 3 



B -frame 



Figure 5: A realization of I and its abstract structure. The pictures do not show an actual example, but 
rather are meant only to illustrate the general structure of a realization (even for our small 6-vertex graph, 
a full matrix would have 28 x 28 blocks). In this particular realization all A-mirrors are perfect mirrors 
associated with a = 101011. This 0-1 vector encodes the vertex set {1, 3, 5, 6}. The three shown edge verifiers 
check that edges (3,5), (4,5) and (1,4) are covered by that vertex set. 

We partition the side-diagonal block h into 

A-mirror h upper left n x n corner 

B-corner b lower right 2x2 corner 

v-screen b remaining entries in column n + 2 

h-screen b remaining entries in row n + 2 

v-gutter b remaining entries in column n-\- 1 

h-gutter b remaining entries in row n + 1. 

Our previous observation implies that the B-corner is completely beige (as it is part of the B-frame), 
and the A-mirror does not contain any beige atoms. The gutters are completely empty and the 
screens completely filled, since 



^b(n+2)+n+l + ''b(n+2)+n+l 
' b{n+2)+n+2 ^ ' 6(n+2)+n+2 



L — n 
L 



„A I B 

*fe{n+2)+n+l ^ ^b{n+2)+n+l 
„A I B 

^b{n+2)+n+2 ^ '^b{n+2)+n+2- 



We define 0-1 vectors P ,cy.^,... ,(3 ,q;™"^^^ to represent positions of the beige atoms in 

the screens: a^ = 1 if and only if i-th atom (from top) in v-screen b is beige and /3f = 1 if and only 
if i-th atom (from left) in h-screen b is beige. 

Since the main-diagonal blocks have no azure atoms, each side-diagonal blocks 6 is a realization 
of ASM[a , P ). Moreover, the main-diagonal blocks and mJ are realizations of, respectively, 
BSM[7v,P ) and BSM[cx."^'^~^,7v), and each other main-diagonal block 5, for b = 1, . . . ,mJ — 1, is 
a realizations of BSM{a ,P ). Lemmas Isl and |6| imply inequalities (H). ■ 
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6.3 The Correctness Proof 

Theorem 1 The problem 3-CCP is NP-complete in the strong sense. 



Proof Clearly, 3-CCP is in NP. To justify the correctness of the reduction described in Section 6.1 
we need to prove that G has a vertex cover of size K if and only if Z is consistent. 

(=^) Suppose that C/ is a vertex cover of size K in G. Define 7 by 7^ = iff u € U . Let T be 
a matrix of the form (^) in which a* = 7, and /3* = 7 for i = 0, . . . .iinaJ — 1. We have 
TT ^ 7 ^ TT . By Lemma Q, T is a realization of I^^. Since U is a vertex cover, Lemma 
implies that T can be extended to a realization of 2. 

(<^) Let T be any realization of I. By Lemma ^, T restricted to azure and beige atoms has the 
form @ . Lemma y implies that the sequence 

(where (3 = tt) has at most K{n — K) + 1 distinct vectors, and thus /3 = (3 for 

some < a < J — 1. By (^), we get 

Define U = {u : a"™" = 0}. Using Lemma 0, we obtain that C/ is a vertex cover. 

To complete the proof, note that Z consists of 6L = O(n^) numbers each bounded by L, so 
the unary encoding of X has size 0{'n}^). Moreover, this encoding can be computed in polynomial 
time. We conclude that 3-CCP is strongly NP-complete. ■ 

7 Final Comments 

We proved that c-CCP is NP-complete for c > 3. Since it is known that 1-CCP can be solved 
efficiently in polynomial time (see [|^), the only unresolved case is for c = 2. 

Relation to multicommodity flows. Consider the following problem: given a bipartite directed 
graph H = (U, V, E), where E is the set of arcs directed from U to V, with each arc having capacity 
1, we want to ship two commodities from the vertices in U to the vertices in V, according to the 
given supplies in U and demands in V. More specifically, for each vertex Uj G C/ we are given a 
supply x^ of commodity a, and for each vertex Vj G ^ we are given a demand y^ of commodity 
a, where a G {1,2}. We wish to compute an integral 2-commodity flow from U to V oi maximum 
total value. Let us call it 2- Commodity Integral 2-Layer Flow, or 2-CI2LF. It is known (see, Q) 
that the 2-commodity integral flow problem is NP-hard for directed networks. We can improve it 
to the 2-layer case. By modifying the argument outlined in Section |^, it is not difficult to show 
that 2-CI2LF is NP-hard as well: simply note that all but two atom types have unique realizations 
which are independent of the given instance G, K of Vertex Cover, and associate the entries not 
occupied by these atoms with the edges of the resulting graph H. (Another proof can be obtained 
by modifying the proof in Q in a similar fashion.) 

The argument above does not imply that 2-CCP is NP-complete, since the graphs corresponding 
to the 2-CCP problem are complete bipartite graphs. This leads to the following open problem: 
Can 2-CI2LF be solved in polynomial time for complete bipartite graphs? 
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Consequences to data security problems. Similar to ^, our result has some consequences 
for problems arising in statistics and data security. 

The reconstruction problem for contingency tables is similar to the 1-CCP problem, except that 
now we allow a realization to contain any non-negative integers (not just O's and I's). Our result 
implies that this problem is NP-hard even when we want to reconstruct a table whose entries are 
in the set {0, 1, /i, //^}, for some given jj.. (To see this, modify the proof by representing each table 
entry in a /i-ary notation, where fi = L + 1, and associate color sums with the coefficients of 1, /x 
and //^.) 

A related problem, arising in the 3D statistical data security problem, is to reconstruct a 3D 
table from its projections, which are called the row, column and file sums. Irving and Jerrum |7| 
proved that this problem is NP-hard even when all file sums are either or 1. The work in |4| 
implies that the problem is NP-hard for L x L x 7 tables and all file sums equal 1. Our result 
improves this result further to tables of size L x L x 4. 
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